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We investigate the Bose-Einstein condensation (BEC, superfluidity) of particle-hole pairs in ul- 
tracold Fermionic atoms with repulsive interactions and arbitrary polarization, which are trapped 
within optical lattices. In the strongly repulsive limit, the dynamics of particle-hole pairs can 
be described by a hard-core Bose-Hubbard model. The insulator - superfluid and charge-density- 
wave(CDW) - superfluid phase transitions can be induced by decreasing and increasing the potential 
depths with controlling the trapping laser intensity, respectively. The parameter and polarization 
dependence of the critical temperatures for the ordered states (BEC and/or CDW) is discussed 
simultaneously. 



c3 
O 



> 
in 

o 
m 
o 

o 

o : 

>>: 
^ : 

Or 



X 



PACS numbers: 03.75.Ss, 32.80.Pj, 71.30.+h 

In recent years, the demonstration of Bose-Einstein 
condensation of particle-particle pairs in homogeneous or 
confined two-component (spin-1/2) ultracold fermionic 
atoms has triggered great theoretical and experimental 
interest. The BCS-BEC crossover in ultracold Fermi 
atomic gases near a Feshbach resonance has been pre- 
dicted by using the resonance superfluidity theory 
and been observed in experiments 2] . For the fermionic 
atoms trapped within optical lattices, the s-wave or d- 
wave particle-particle pairs can undergo a phase transi- 
tion to a superfluid state when the inter-component in- 
teraction is attractive or repulsive 

With the mechanism of the superfluidity of atom-atom 
pairs in ultracold Fermi atomic gases being explored more 
and more deeply, the question arises whether the atom- 
hole pairs in ultracold Fermi atomic gases can undergo 
a BEC phase transition similar to electron-hole pairs 
0. Firstly, due to the essentially non-equilibrium na- 
ture of the system of particle-hole pairs, which usually 
does not appear in the system of particle-particle pairs, 
the particle-hole system becomes an ideal system for ex- 
ploring the non-equilibrium quantum mechanics at the 
frontier of many-body physics. Additionally, since all 
condensed particle-hole pairs can emit photons in tan- 
dem, the quantum coherence in such a condensate will 
reveal novel optical effects and nonlinear optical dynam- 
ics, which can not be shown by the condensate of particle- 
particle pairs. This provides possible applications in ul- 
trafast digital logical elements and quantum computa- 
tion. Furthermore, the controllable interaction strength 
and kinetic energy (or hopping strength) in the atomic 
systems open up the very exciting potential to investigate 
macroscopic quantum coherence and superfluidity under 
some extreme conditions that never exist in traditional 
electronic systems. 

In this letter, we show that the atom-hole pairs in ar- 
bitrarily polarized spin-1/2 ultracold Fermi atoms with 
repulsive interaction, which are trapped within optical 



lattices, can undergo a superfluid phase transition sim- 
ilar to the ultracold bosonic atoms confined in optical 
lattices In the strongly repulsive limit, the dynamics 
of atom-hole pairs can be described by a hard-core Bose- 
Hubbard model. Then, the phase transition is analyzed 
with the derived Bose-Hubbard model. At the same time, 
the critical temperature for the ordered states (charge- 
density- wave and/or Bose-Einstein condensation) is dis- 
cussed within the mean-field theory. 

Consider an ensemble of arbitrarily polarized ultracold 
fermionic atoms occupying two different hyperfine states 
15*) and |P), which are trapped within the optical lattices. 
For simplicity, we assume the optical lattice potentials 
as Vq' p (x) — J2j=i Vq' p cos 2 (kxj). The wave vector k is 
determined by the laser wave-lengths, d (=1, 2 or 3) is 
the dimension of the optical lattices, and V^ ,p are pro- 
portional to the laser intensity. For sufficient low tem- 
perature, all atoms will be localized into the lowest Bloch 
band, and the system can be described by an asymmetric 
Fermi-Hubbard Hamiltonian [3j 



H = - J2 (tsfsifsj +tpfpifpj) 
(i,3) 

+ J2( e sn si + e p n pi ) + UJ2 n si n pi . 



(1) 



Here (fai) are fermionic creation (annihilation) oper- 
ators for localized atoms in state \a) on site i, n a i = 
f a i. The symbol represents summing over the 

nearest-neighbors and e s (e p ) is the single-atom energy 
of the atoms in state \S) (\P)). The state-dependent 
tunneling (hopping) t s (t p ) between nearest neighbors can 
be induced by varying the potential depth Vq(Vq) with 
controlling the laser intensity 0. Usually, the tunneling 
strengths increase with the decrease of potential depths. 
Below the unitary limit, the on-site interaction U is pro- 
portional to the s-wave scattering length between atoms 
occupying different hyperfine states. The s-wave scatter- 
ing between atoms occupying the same hyperfine state is 
absent due to the Pauli blocking. 
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The average number of atoms per site (filling number) 
n and the polarization 7 of the considered system are 
defined as 



(n si + n pi )/N L 



7 



)/ }Xn S j + n P i). 



(2) 



(3) 



The symbol JVl is the total number of lattice sites. In 
the following, we focus our interests on the half-filled case 
(n = 1), i.e., one atom per site. 

The ground state energy per atom depends upon both 
the polarization and the energy difference (Ae = e p — e s ) 
between two occupied states. With the definition of po- 
larization, the ground states can be divided into five dif- 
ferent regimes: non-polarized (NP) ground states with 
7 = 0, partially polarized in state \S) (PPS) with 
< 7 < 1, partially polarized in state \P) (PPP) with 
-1 < 7 < 0, fully polarized in state \S) (FPS) with 7=1 
and fully polarized in state \P) (FPP) with 7 = — 1. 
For the one dimensional lattices (d — 1) with state- 
independent hopping (t s = t p = t), these regimes can 
be exactly obtained with the Bethe-ansatz 0. The Ae 
has two critical values 



and 



M-2*. 



Ji (w)dw 



lu[ l +eX p(l^)] 



for U < 0, 
for U > 0, 




u 
2 



for 
for 



U < 0, 

u> 0, 



(4) 



(5) 



corresponding to the boundaries between different 
regimes. Here, Ji(w) is the first kind Bessel function 
with first order. The non-polarized, partially polar- 
ized and fully polarized regimes satisfy |Ae| < Ae^, 



|Ae < 



Ae§ and 



|Ae| > AeJS respectively (see 



Aef < 
Fig. 1). 

In the strongly repulsive limit (0 < t a>p <C U), the 
Fermi-Hubbard model is equivalent to an effective spin- 
1/2 Heisenberg model |8j. For the case of infinitely re- 
pulsive limit (U/t s>p — > +00), the ground states (lowest 
energy states) have only one atom for each site, and their 
charge degrees of freedom are frozen. Under the strongly 
repulsive condition, U/t StP ^> 1, one can introduce the 
bosonic operators 



b~j & f s jfpj, bj ^> /+/ s 
rij = b+bj \ + i(n aj - 



(6) 



°P3 > 



for the atom-hole pairs on site j. The operator bj(bj) 
creates (annihilates) a pair of S-atom (atom in 15*}) and 
P-hole (hole in |P)) on site j. These operators with differ- 
ent lattice indices are commutable. However, to exclude 
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FIG. 1: Left: Ultracold fermionic atoms in one-dimensional 
optical lattices with half-filling and state-independent hop- 
ping. The dots and circles denote the atoms and holes 
(no atoms) respectively. Right: Polarization regimes of the 
ground states for the one-dimensional lattices with state- 
independent hopping. 



the multiple occupation at each lattice which comes from 
Pauli blocking, the operators with same lattice indices 
have a property like Fermi particles. In other words, 
the interaction between bosons on the same lattice is in- 
finitely repulsive. Using the perturbation theory devel- 
oped by Takahashi y| (b and b + correspond to er~ and 
cr + ), up to third order terms of the perturbation param- 
eters (hopping strengths), we obtain the atom- hole pairs 
obey the hard-core Bose-Hubbard Hamiltonian 



(i,3) 



(i,3) 



(7) 



Denoting t p = od s — at, we obtain the hopping strength 
J = 4tpt s /U = 4at 2 /U, the nearest neighbor interaction 
strength V = 2(t 2 a +t 2 p )/U = 2{l + a 2 )t 2 /U and the chem- 
ical potential fi = Ae + ZV/2 = Ae + Z(l + a 2 )t 2 /U . For 
the cubic lattices, the total number of the nearest neigh- 
bors Z equals 2d. The above hard-core Bose-Hubbard 
model can be mapped onto an anisotropic spin-1/2 XXZ 
Heisenberg model with = J , J z = V and an effective 
magnetic field B z = Ae |2(. The antiferromagnetic-Z or- 
der, XY-order and fully magnetized states in XXZ model 
correspond to the CDW phase, superfluid phase and fully 
polarized insulator phase of the atom-hole pairs, respec- 
tively 0, . The superfluidity means the Bose conden- 
sation of atom-hole pairs in their momentum spaces. 

At zero temperature, the ground states for the atom- 
hole pairs have three different phases: (i) charge-density- 
wave phase similar to a solid phase with zero polar- 
ization (7 = 0) corresponds to the half-filled case of 
the hard-core Bose-Hubbard model ((b + b) = 1/2), (ii) 
Bose-Einstein condensation phase with non-zero super- 
fluid order parameter (6), and (iii) insulator phase with 
the largest polarization (I7I = 1) corresponds to the 
empty ((b+b) = 0) or the fully-filled ((b+b) = 1) case 
of the hard-core Bose-Hubbard model. The difference 
between superfluid and insulator phases indicates that 
it need a non-fully polarized atomic gases to support 
the atom-hole BEC. From the equivalence between the 
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hard-core Bose-Hubbard model and the spin- 1/2 XXZ 
Heisenberg model, using the path- integral method 
one can obtain that the fully polarized insulator phase 
appears when |Ae| /U > Z{ t/U) 2 (l + a) 2 , the BEC 
phase exists if (Z/2)(t/U) 2 y/(l - a 2 ) 2 < |Ae| /U < 
Z(t/U) 2 (\ + a) 2 , and the CDW phase emerges when 
\Ae\/U < (Z/2)(t/U) 2 y / (l-a 2 ) 2 . The separatrix be- 
tween CDW phase and BEC phase corresponds to a 
first order phase transition. The points on this sepa- 
ratrix means the coexistence of both phases, they repre- 
sents the supersolid phase. These conditions also show 
the CDW - superfluid and insulator - superfluid tran- 
sitions occur at |Ae| /U = {Z /2){t/U) 2 ^/(l - a 2 ) 2 and 
|Ae| /U = Z{t/U) 2 {\ + a) 2 , respectively. 

In FIG. 2, we show the phase diagram for lattices 
of arbitrary dimensionality with hopping ratio a = 2. 
For fixed values of hopping ratio a, energy difference 
Ae and on-site repulsive interaction strength U, increas- 
ing (decreasing) the hopping strength t will induce an 
insulator-superfluid (solid-superfluid) transition. This 
means that Bose condensation of the atom-hole pairs ex- 
ists for mediate hopping strengths. For larger (smaller) 
hopping strength, the ground states fall into the phase 
of CDW (fully polarized insulator) . In the case of state- 
independent hopping, a — 1, the CDW region becomes 
a line localized at Ae = 0. This is consistent with the 
results of the antiferromagnetic phase in Refs. Q. 
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FIG. 2: Zero-temperature phase diagram of the ground states 
for the atom-holes in arbitrarily dimensional lattices with a = 
2. 

At finite temperatures, due to the thermal fluctuations, 
the ordered phases will be destroyed when the temper- 
ature is above some critical temperatures. Within the 
framework of the mean-field theory [l2T |. a continuous 
phase transition between the CDW and the normal liq- 
uid (NL) takes place at the critical temperature 



rpC 

1 CDW 



Z 



(l+a 2 )t 2 
U 



(1-7 2 ), 



(8) 



and a similar phase transition between the superfluid and 
the normal liquid occurs at 



1 SF — ~ — 



Z 2at 2 



7 



U arctanh(j) 



(9) 



Here, ks is the Boltzmann constant. The bicritical 
polarizations 7 = ±-f BC (a)( r y BC > 0) are given by 
Tg F = Tq DW and | 7 |^ 1. Below the critical tempera- 
tures, there are two coexistence regions of CDW and su- 
perfluid, which correspond to the supersolid regions. The 
boundaries between the superfluid and supersolid and be- 
tween CDW and supersolid can be obtained by using the 
Landau expansion At zero temperature, the critical 
polarization corresponding to the superfluid-supersolid 
transition is given as ^ SFC = |(1 — a)/(l + a)\. For 
the case of state- independent hopping (a = 1), the criti- 
cal polarization -f SFC = 0, it means that the CDW/solid 
and supersolid regions shrink to a line localized at 7 = 0. 
This indicates that the atom-hole BEC in non-polarized 
atoms with state-independent hopping has the highest 
critical temperature for the superfluid phase. 

The finite temperature phase transitions rely on the 
hopping ratio a and the polarization 7. For state- 
independent hopping (a = 1), CDW - NL and super- 
fluid - NL transitions occur in non-polarized (7 = 0) and 
polarized case (7 / 0), respectively. For state-dependent 
hopping (a 7^ 1), the transition routes become more com- 
plex. The CDW - NL, supersolid - CDW - NL, super- 
fluid - supersolid - CDW - NL, and superfluid -NL phase 
transitions take place when 7 = 0, < I7I < j SFC , 
rySFC < |/y| < 7 BC , and 7I > 7 SC , respectively. 

The critical temperatures are determined by both the 
parameters and the polarization. The parameter depen- 
dence is similar to the one of Refs. 3], T c cx t 2 /U. 
Thus, to increase the critical temperatures, one has to de- 
crease the potential depths Vq' p to increase the hopping 
strengths. The polarization dependence of the critical 
temperatures and the finite temperature phase diagram 
are shown in FIG. 3 for the state-dependent hopping case 
with a — 2. 

The previous consideration includes only the terms up 
to the third order of the perturbation parameter. Includ- 
ing up to the fifth order terms, the Hamiltonian (7) reads 
as 



H B = -/iE"i + JiE bfbj + ,h E b+h 

i (id) ((hk)) 

+Vl J2 n i n i +V2 J2 n i n k- 
(id) ((i,k)) 



(10) 



Here, ((i, k)) represents summing over the next-nearest- 
neighbors. The parameters are determined by /i = Ae + 



4t s t„ 



u 3 



and Vo — 



},Vi 



Z{V 1 + V 2 )/2, Ji 

2U 2 t s t P 2 

4tgtp)/(2t 2 tp). The corresponding critical temperatures 
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FIG. 3: Mean-field finite-temperature phase diagram for the 
atom-hole pairs in state-dependent hopping case with a = 2. 

are formulated as 

r&w = ^ • (Vx - v 2 ) ■ (i - 7 2 ), (ii) 

T * F = i ' (Jl + J2) • arrfdhW (12) 

Because Ji <C J\ « J and V2 <C Vi » V", the above 
equations indicate that the critical temperatures are only 
shifted a little bit by the high-order terms. 

In summary, we have demonstrated the existence of 
Bose-Einstein condensation of atom-hole pairs in arbi- 
trarily polarized ultracold fermionic atoms confined in 
optical lattices with half filling. In the strongly repul- 
sive limit, the atom-hole pairs obey a hard-core Bose- 
Hubbard Hamiltonian. For a polarized insulator phase, 
the particle-hole pairs undergo an insulator - superfluid 
transition when the hopping between nearest neighbors 
is increased. For a charge-density-wave (CDW) phase, 
the pairs undergo a CDW - superfluid transition when 
the hopping is decreased. The mean-field results indi- 
cate that the finite temperature phase transition depends 
upon not only the system parameters but also the polar- 
ization. 

To realize the considered model, one can prepare an 
ultracold two-component atomic Fermi gases with arbi- 
trary polarization |2| , then load them into an optical lat- 
tice with one atom per site. The optical lattices can be 
produced with a series of standing-wave lasers. Similar to 
the experimental realization of Tonks gas [l4j (however 
the physical details are only loosely related) , the strongly 
repulsive limit U/t SiP 3> I can be reached by increasing 
the s-wave scattering length with Feshbach resonances 
0| (but still below the unitary limit 0) and/or by de- 
creasing the hopping strengths t S)P with controlling the 
laser intensity. The applied magnetic field will also in- 
duce a energy difference between two occupied levels due 



to the Zeeman effects. To observe the superfluidity of 
the atom- hole pairs, one can use Bragg scattering ap- 
proach to detect the elementary excitations spectrum of 
cold atoms by monitoring the scattered atoms versus the 
frequency difference between two lasers |l7| . which form 
the light grating. In the superfluid phase, it will appear a 
distinguished peak corresponding to the collective Bigoli- 
ubov quasi-particle excitations. In a condensed system of 
interacting atom-hole pairs within an applied electromag- 
netic field, the stimulated two-photon emission (second 
order harmonics emission) process can also give a credi- 
ble evidence for the atom- hole BEC . In this process, 
the number of the condensed atom-hole pairs decreases 
with the emission of the second order harmonics of the 
applied light. 

The author C. Lee acknowledges useful discussions 
with professor Peter Fulde and Dr. Joachim Brand. This 
work is supported by foundations of Max Planck Insti- 
tute for the Physics of Complex Systems. 



* Electronic addresses: chlee@mpipks-dresden.mpg.de; 
chleecn@hotmail.com 

[1] M. Holland, et al., Phys. Rev. Lett. 87, f20406 (2001); M. 
L. Chiofalo, et al., Phys. Rev. Lett. 88, 090402 (2002); 
Y. Ohashi and A. Griffin, Phys. Rev. Lett. 89, f30402 
(2002); L. D. Ca rr, et al., |cond-ma t/0308306 (2003); J. 
Kinnunen, et al., ||cond-mat /040f 543 (2004^ 

[2] S. Jochim, et al., Science 302, 2101 (2003); C. A. Regal, 
M. Greiner and D. S. Jin, P hys. Rev. Lett. 92 040403 
(2004); M. Bartenstein, et al., |cond-m at/0401109 (2004). 

[3] W. Hofstetter, et al., Phys. Rev. Lett. 89, 220407 (2002); 
M. Rigol, et al., Phys. Rev. Lett. 91, 130403 (2003). 

[4] L. V. Keldysh and A. N. Kozlov, Soviet Phys. JETP 27, 
521 (1968); D. S. Chemla and J. Shah, Nature 411, 549 
(2001); I. E. Perakis, Nature 417, 33 (2002); L. V. Butov, 
et al., Nature 417, 47 (2002); C. W. Lai, et al., Science 
303, 503 (2004). 

[5] D. Jaksch, et al., Phys. Rev. Lett. 81, 3108 (1998); H. T. 
Stoof, Nature 415, 25 (2002); M. Greiner, et al., Nature 
415, 39 (2002), A. Smerzi, et al., Phys. Rev. Lett. 89, 
170402 (2002), C. Menotti, A. Smerzi, and A. Trombet- 
toni, New J. Phys. 5, 112 (2003); B. Damski, et al., Phys. 
Rev. Lett. 91, 080403 (2003). 

[6] E. Altman, W. Hofstetter, E. Demler and M. D. Lukin, 
New Journal of Physics 5, 113 (2003). 

[7] E. H. Lieb and F. Y. Wu, Phys. Rev. Lett. 20, 1445 

(1968) ; M. Takahashi, Prog. Theor. Phys. 42, 1098 

(1969) ; C. Yang, A. N. Kocharian and Y. L. Chiang, 
J. Phys.: Condens. Matter 12, 7433 (2000). 

[8] M. Takahashi, J. Phys. C: Solid State Phys. 10, 1289 

(1977); G. Fath, Z. Domahski, and R. Lemahski, Phys. 

Rev. B 52, 13910 (1995); C. D. Batista, Phys. Rev. Lett. 

89, 166403 (2002). 
[9] T. Matsubara and H. Matsuda, Prog. Theor. Phys. 16, 

569 (1956); G. G. Batrouni, et al, Phys. Rev. Lett. 74, 

2527 (1995); G. G. Batrouni and R. T. Scallettar, Phys. 

Rev. Lett. 84, 1599 (2000); G. Schmid, et al., Phys. Rev. 

Lett. 88, 167208 (2002). 



5 



[10] L. J. de Jongh, Physica C 161, 631 (1989). 

[11] S. Sachdev, Quantum Phase Transitions (Cambridge 

University Press, Cambridge, England, 1999), Chap. 13; 

F. Hebert et al., Phys. Rev. B 65, 014513 (2002); A. Auer- 

bach, F. Berruto and L. Capriotti, cond-mat/9801294 

(1998). 

[12] S. Robaszkiewicz, R. Micnas and K. A. Chao, Phys. Rev. 

B 23, 1447 (1981); A. Aharony ans A. Auerbach, Phys. 

Rev. Lett. 70, 1874 (1993). 
[13] K. S. Liu and E. Fisher, J. Low Temp. Phys. 10, 655 

(1973). 

[14] B. Paredes, et al., Nature 429, 277 (2004). 
[15] E. Timmermans, et al., Phys. Rep. 315, 199 (1999); S. 
Inouye, et al., Nature 392, 151 (1998); J. L. Roberts, 



et al., Phys. Rev. Lett. 85, 728 (2000); K. Dieckmann, 
et al, Phys. Rev. Lett. 89, 203201 (2002); T. Loftus, 
C. A. Regal, C. Ticknor, J. L. Bohn, and D. S. Jin, 
cond-mat/01 11571 (2001). 
[16] S. J. J. M. F. Kokkelmans, et al., Phys. Rev. A 65, 053617 

(2002) ; R. Combescot, Phys. Rev. Lett. 91, 120401 

(2003) ; S. Gupta, et al., Phys. Rev. Lett. 92, 100401 

(2004) ; T. -L. Ho, Phys. Rev. Lett. 92, 090402 (2004). 
[17] J. Stenger, et al., Phys. Rev. Lett. 82, 4569 (1999). 

[18] Yu. E. Lozovik and A. V. Pushnov, JETP 88, 747 (1999); 
T. Portengen, Th. Ostreich, and L. J. Sham, Phys. Rev. 
Lett. 76, 3384 (1996); T. Portengen, Th. Ostreich, and 
L. J. Sham, Phys. Rev. B 54, 17452 (1996). 



